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Abstract 

Texture zeros in the quark Yukawa matrices generally lead to precise and simple ex- 
pressions for CKM matrix elements in terms of ratios of quark masses. Using the new 
data on b— decays we test a particularly promising texture zero solution and show that 
it is at best approximate. We analyse the approximate texture zero structure and show 
it is consistent with experiment. We investigate the implications for the CKM unitarity 
triangle, measurements at BaBar and BELLE as well as for the theories which invoke 
family symmetries. 



1 Introduction 

The structure of quark and lepton mass matrices provides us with a rare insight into the physics 
beyond the Standard Model which may directly probe the underlying theory at the gauge 
unification or Planck scale. While the quark mass matrices and the CKM matrix, V , 
are intimately related, measurement of the eigenvalues of the mass matrices and the matrix 
elements of V is not sufficient to determine the structure of the full mass matrix and 
of the matrix of Yukawa couplings giving rise to them. Given this under-determination, the 
phenomenological approach most often used is to make some assumption about this structure 
and explore the experimental consequences for the V$ . A particularly promising starting 
point assumes that there are anomalously small entries in the up and down quark Yukawa 
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matrices - "texture zeros''^. These lead to relations for the in terms of ratios of quark 

masses which do not involve any unknown couplings and hence can be precisely tested. Various 
texture zeros have been studied. For the case of symmetric mass matrices a systematic analyses 
determining which combinations of textures involving 4, 5 or 6 zeros for the U, D matrices are 
compatible with data was carried out in Q. The main reason for looking for such texture zero 
solutions is that they may shed light on physics beyond the Standard Model, for example the 
presence of a new family symmetry relating different generations. 

The new generation of b-factory experiments has led to more precise measurements of the 
CKM matrix elements that, together with the progress in understanding hadronic uncertain- 
ties [0 and light quark masses, allows us to test texture zero structures to a greater precision 
than has hitherto been possible. In this paper we will study the most promising structure 
based on simultaneous zeros in the U and D mass matrices at the (1,1) and (1,3) positions. In 
the additional hypothesis of equal magnitude of the (1,2) and (2,1) entries and of sufficiently 
small (3,2) entry, this structure gives three precise relations between and ratios of quark 

masses, leaving only one CKM element undetermined. We find that the new experimental and 
theoretical information suggests that (at least) one of the hypothesis leading to those precise 
texture zero relations needs to be relaxed. The simplest possibility is that while the (1,3) ele- 
ment is small it is non- zero so that the texture zero is only approximate. As a result one of 
the relations, the texture zero prediction for | V^&/V^b| , is modified, as is suggested by the new 
precise data; the other two relations are less affected. 

In the context of an underlying family symmetry this result is to be expected for the family 
symmetry usually requires texture zeros to be only approximate and, in some cases, actually 
predicts the order at which the approximate zero should be filled in. For example, a very simple 
Abelian family symmetry predicts the (1,3) element should be nonzero at a level consistent with 
the new data while requiring the (1,1) element should be much smaller, preserving the remaining 
two texture zero predictions. 

Another possibility is that the (3,2) entry in the down quark mass matrix is not as small 
as the (2,3) entry (barring cancellations, the latter must be smaller than the (3,3) entry by a 
factor O {\V c b\)). Such an asymmetry can also be easily achieved in the context of Abelian or 
non- Abelian family symmetries and offers an intriguing connection with neutrino physics. A 
sizeable (3,2) entry in the down quark mass matrix is in fact a generic prediction of a class 
of unified models of quark and lepton masses and mixings. In such models, a large leptonic 
mixing angle accounting for the atmospheric neutrino anomaly originates from a sizeable (2,3) 
entry in the charged lepton mass matrix. The GUT symmetry then forces the (3,2) element in 
the down quark matrix to be of the same order of magnitude [|]]. However we note that the 
original symmetric structure can also lead, very naturally, to a large neutrino mixing angle ||. 
Finally, the third situation leading to a correction to the texture zero relations is that the entries 

1 Strictly texture zeros can only apply at a single mass scale (the GUT or string scale?) and will be filled 
in by Renormalisation Group running. However, in general, such effects are very small and the texture zeros 
persist to a good approximation at all scales. 
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(1.2) and (2,1) in the up quark mass matrix are not equal in magnitude. This can happen in 
non unified abelian models due to different order one coefficients. We will not investigate this 
possibility in this paper because it destroys one of the successful texture zero predictions (see 
the discussion in Section ^ below). From this one sees that detailed tests of the texture zero 
relations will help to identify the underlying family symmetry. 

Our analysis does not take into account possible new physics contributions to the pro- 
cesses constraining the CKM parameters. Such contributions might affect the experimental 
determination of |V^/14 S | through modifications to B mixing and the CP- violating part of K 
mixing (or both) and consequently affect the corresponding texture zero relation. However 
we do not expect IK^/V^I to be similarly affected where the main uncertainty instead lies in 
the hadronic modelling of charmless semileptonic B decays. Our analysis addresses the phe- 
nomeno logical problem of the jV^/V^I prediction of texture zero structures and its remedies 
and, as a consequence, most of our analysis would not be affected by new physics effects. Very 
recent measurements of sin 2(3 from BaBar and BELLE || do raise the possibility of beyond 
Standard Model contributions to CP- violation and therefore the quantitative fits performed in 
Sections § and Q might be affected by supersymmetric contributions to K and B mixing. We 
briefly discuss how low values of sin 2(3 affect our analysis. 

The paper is structured as follows. In Section |2| we present the experimental tests of the 
texture zero predictions following from zeros in the U and D mass matrices at the (1,1) and 

(1.3) positions. In Section |3] we discuss the implications for the mass matrices following from 
the need to modify one of the texture zero relations. We develop a perturbative expansion 
which allows us to identify the possible corrections to the texture zero predictions. In Section 
|] we consider the implications of such structure for an underlying family symmetry and finally 
in Section 5 we present conclusions. 

2 Tests of texture zero predictions 

In what follows we assume that the off-diagonal entries are small relative to their on-diagonal 
partners so that one may develop a perturbative expansion for the CKM matrix elements 
j7|, ||, [| fTO] . This is a reasonable starting point because it immediately leads to small mixing 
angles consistent with observation. We further assume here that there are texture zeros in the 
(1,1) and (1,3) elements, that the (1,2) and (2,1) elements have equal magnitude and that the 
texture is approximately symmetric ( (3,2) ~ (2,3)). These assumptions lead to the texture zero 
relations [|, |, 0, 0, 
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Fixed Parameters 



Parameter 
Gp 
M w 

fx 
m K 
ArriK 
\£k\ 

V2 

m Bd 

Vb 



Value 

1.16639 x l(T 5 GeV- 2 
(80.42 ± 0.06)GeV 
(0.161 ±0.0015)GeV 
(0.497672 ±0.000031)GeV 
(3.491 ± 0.009) x 10 15 GeV 
(2.271 ±0.017) x 10~ 3 
(0.574 ±0.004) 
(5.2792 ±0.0018)GeV 
0.55 ±0.01 

(5.3693 ±0.0020)GeV 



Reference 







m 







El 



IH 



m 

pi 











Table 1: Fixed Parameters. 



We will prove that, quite generally, is approximately the Standard Model CP violating phase 
(for more restricted cases see 0, |14|1). 

Barbieri et al (BHR) |]TTJ emphasized that these relations lead to a very tight determination 
of the CKM unitarity triangle. In terms of the re-scaled Wolfenstein parameters, p = cp, 
fj = ct], c = (1 - A 2 /2), the ratios IK&I/IK&I and |T4 rf |/|V^ s | are 



IK 



ub 



W, 



cb\ 



p2 _|_ ^2 



\V t 



td\ 



ts 



1 - p) 2 + r] 



n'2 



(3) 



In the context of the Standard Model, the measurable quantities which give information on 
p and f\ are 

(a) the ratio |V^&|/|V^,| obtained from semi-leptonic decays of B mesons, 

(b) Am^ and Amg s which are the mass differences in the B® — B° d , B® — B° s systems, 

(c) \sk\ the parameter related to CP violation in the K, K system, and 

(d) sin 2(3 (/3 is one of the angles in the unitarity triangle) obtained from CP asymmetries 
in various B decays. 



2.1 Standard Model (SM) fit 

In comparing the texture zeros with experiment we proceed in two stages. We first use the 
latest data to find A, p and fj and then compare the result with eq@. Our procedure is 
to construct a two-dimensional probability density for p and fj [|17| from the constraints of the 
above measurements. Entering in the fits are several parameters which have been well measured 
and which we choose not to vary. These are given in Table [TJ. At present we have only an upper 
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limit for Am^ and we employ the so-called 'amplitude method' to include this information 



into the fit[H8 



Our fit assumes the Standard Model (SM) relations between the experimental measurables 
and the CKM matrix elements. In a supersymmetric extension of the Standard Model we 
expect that there will be corrections to these measurables; |e#| is particularly sensitive to such 
corrections. To allow for this possibility we carry out a separate fit in which the data on \ejc\, 
together with that on sin 2/3 for which there are rather different measurements, is dropped. The 
mass differences ArriB d , Ams s and sin 2/3 might also be affected by supersymmetry, especially 
if the structure of squark mass matrices is determined by the flavour symmetry accounting for 
quark masses and mixings. We do not consider this possibility here. 

The formulas used in the standard model fit are 



Am Bd = C A A 2 \ b [(l - pY + ^K4B Bi ^S( 



x 



(4) 



6tt 2 



-. Here S(x*) is the standard Inami-Lim function |19| and f B B B is the 



where C AmB 

product of the B meson decay constant and the B parameter analogous to Bk in the K system 

.2 1 



where 



Am Bs = Amn 



< 



m Bd A 2 (1 - p) 2 + ?7 2 ' 



[b s \/B Bs 
ft 



B B d 



(5) 



(6) 



For the 6k parameter we have 

: C e B K A 2 X 6 fj [-^x* + A 2 \ A (l - p - (p 2 + f - p) A 2 ) ti*3{x* 



\£k\ 



(7) 



G F^K mKm w _ The short distance QCD corrections are contained in the coefficients 



where C e 

7]*, which have been computed at the Next to Leading-logarithmic Order (NLO) in p0 |. The 



NLO calculation requires the use of the one-loop relation between the pole mass and the run- 
ning mass mf oZe = m\ (l + ""^ ^ f) i n the MS. The coefficients rj* have been evaluated at 
A*ff=37lMeV. All the starred quantities in Table | are given in terms of these quantities, 
for example x* = [m*] 2 /M^. 

Finally the angles of the unitarity triangle are given by 



sin 2/3 



sin 2a 



sin 27 



27/(1 -p) 
f] 2 + (1 - p) 2 

2fj (r/ 2 + pjp - 1)) 

in 2 + (1 - p) 2 ) {n 2 + p 2 ) 

p2 _|_ fj2 



(8) 
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Fitted Parameters 


1— ' Q y Q TYl nf OT 


Value 


( Q llCOl QTl H QT DTVATO 

VjrcLU.bblcLIl _ -r IcLL tjllUIo. 


Refer en. 


4 


n s^a -+- n n^fi 






* 




\ 

/A 


D 2196 ± D 0023 






* 




It/ iCLEO 
| v ub\ 


32 K x ID -4 


(_lo q i tr r\ y i n-4 


[ 




1 


W J LEP 


41.3 x 10~ 4 


f+fi 3 + 311x1 0~ 4 




22] 


IK. J 

| " CO | 


(41 ± 1 6) x 10~ 3 






1 






0.87 ± 0.143 


0.06 ± 0.13 
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c 


(1.3 ±0.1)GeV 






* 




nit 
t 


(167±5)GeV 






* 




a 


1.38 ± 0.53 








0.47 ±0.04 




[ie] 


Am B d 


(0.487 ±0.014)p S - 1 




[ 




I 


fB d \J B Bd 


(0.230 ±0.032)GeV 


±0.025-±0.020 




§ 




t 


(1.14 ±0.064) 


0.04 ±0.05 




i 




Am Bs > 


ISps- 1 at 95% C.L. 






23] 


sin 2f3 


0.41 ±0.17 










sin 2/3 


0.47 ±0.16 






24 





Table 2: Fitted Parameters. The parameters marked with * have been computed here with the new 
data from [15]. 



Using these expressions we carried out fits to the parameters listed in Table || keeping the 
well determined parameters listed in Table |l| fixed. 

The results in the p — fj plane of the fit are shown in Fig.|l]. The confidence limits shown 
correspond to 68%, 95% and 99%. In the experimental fits (which assume the SM) the results 
of \V u b/V c b\ for both CLEO and LEP collaborations (see Table |2|) were included. Although the 
result from the CLEO collaboration is lower than the result from LEP, they are consistent 
within one a. The combination of both, assuming a Gaussian distribution for the experimental 
errors and a flat distribution for the theoretical errors, gives a value of 0.087 ± 0.010, which is 



consistent with the PDG value of 0.090 ± 0.025 [15 



As can be seen from Fig [I] the constraints on p and fj would be considerably strengthened by 
a measurement of Am^, which presently has only a lower limit, and by an improvement in the 
precision on sin 2/3. We can see from Fig.[T] this agrees within 2a with the the new value for the 
parameter sin 2/3=0.47 ± 0.16 (combining BaBar and Belle results with those from CDF and 
24J). However there is a deviation at the la level which may be a hint for physics 



ALEPH 



beyond the SM. 
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Figure 1: The SM fit of Section 2.1 to (Kb/V^l, Am# s (lower limit), Am^, \ck\ and the recent 
result for sin 2/3. The lines indicate the region of la and, in the case of sin 2/3, also the 2a region, 
demonstrating that the new value is still consistent with the rest of SM constraints within 2a. The 
CL are at 99%, 95% and 68%. 



2.2 Comparison with the texture zero predictions. 

We are now able to compare the experimental results with the texture zero predictions of eq(|I]). 
In Fig.|2| we show the region in the p — fj plane allowed by these relations together with the 
various constraints following from the processes (a)-(d). In this fit we have taken symmetric 
forms for the U and D Yukawa matrices with texture zeros in the (1,1), (1,3) and (3,1) positions. 
Comparison with Fig [I] show that the predictions are hard to reconcile with the data, being 
consistent only at greater than the 99% CL (for the case \ck\ and sin 2/3 are not included). Fig.|| 
shows that one of the reasons for the poor agreement is the measurement of . Given the 
fact the CLEO and LEP measurements differ considerably it is of interest to consider whether 
the discrepancy disappears if we use only the value for the CLEO collaboration. In fact we 
find that this only marginally changes things. One may see that the significant improvement 
in the experimental measurements, particularly of and the improved lower limit on Am# s , 
strongly disfavour this promising texture zero scheme. As remarked before, this conclusion 
holds in extensions of the SM too, since new physics effects might affect the mass difference 
ArriB a , but cannot alleviate the disagreement with V u b/V c b. 

Given this discrepancy, do we have to abandon the texture zero solution completely? In fact 



we do not, as we now show. The problematic relation |\4b/K&| = ym«K (and, to a lesser 
extent, the other two texture zero predictions) depends on three assumptions. 

• Texture zeros: the matrix elements Y\%, Y31, Y\\ are negligibly small both in the up 
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Figure 2: The fit for the relations |Vtd|/|Vts| = y/md/m s and |V^&|/Vd,| = \Jm u /m c . The experi- 
mental constraints + SM interpretation indicate that these relations are disfavoured. 

(Y = U) and down (Y = D) sector. Actually, as mentioned above, the underlying theory 
generating the texture zero is unlikely to guarantee a particular mass matrix element is 
absolutely zero. Moreover an exact zero can only apply at a single mass scale for radiative 
effects necessarily generate contributions to all the elements of the quark mass matrices. 

• Small higher order corrections in the perturbative diagonalization of the up and down 
quark mass matrices. The expected correction is of order 7% or less if \D 32 \ ~ | V^j, 1 1 Z? 3 3 1 
and becomes important for larger values of \D 32 \ ■ Larger values of I-D32 1 can arise in unified 
models in which a large leptonic mixing originates from the charged lepton sector. They 
also arise in models with a texture zero in the (2,2) position and no cancellation between 
down and up quark contributions to Va,. In this case one has in fact ID23/-D33I ~ \Vcb\- 
As a consequence, an asymmetry |Z) 32 / -D33I > I-D23/-D33I is required in order to account 
for the value of m s /m b = \D 23 /D 33 ■ D 32 /D 33 \. 

• \Uu\ = \U 2 i\- This condition is usually met in unified models. In the case of SU(5) models 
for example, Ui 2 and U 2 i are both generated by operators that can be written in the form 
(4>) TiT 2 H, where T 12 are the tenplets of the first and second family containing the up 
quarks, H is the up Higgs fiveplet and represents a (normalized) set of fields whose vev 
is SM invariant. The relation \U\ 2 \ = \U 2 i\ follows unless (4>) breaks SU(5). In the latter 
case, SU(5) Clebsh coefficients will differentiate \Ui 2 \ and \U 2 %\, but often in a too violent 
way, leading to an even worse disagreement with the Kt&/K* prediction. Moreover, the 
analogous operator in the down quark and charged lepton sector would spoil the successful 
relation mdm s /ml ~ m e m M /m^. The condition \Ui 2 \ = \U 2 \ \ is automatically met in some 
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non-Abelian models |2^], |26[. In any case we regard the phenomenological success of the 
relation for \V US \ given by eq(|2|) as a result that should be preserved. 

In what follows we focus on the possibility that either the first or the second assumption is not 
fulfilled. In the first case we still assume a symmetric structure for U and D and since in the 
second case D is manifestly asymmetric, we refer to the these two scenarios as the symmetric 
and asymmetric texture cases. 

In particular, in Section |3.2| we study in detail textures with small but non negligible (1,3) 
element. As we discuss there, the order at which this element arises is a characteristic prediction 
of a family symmetry so determination of this element is a discriminator between various 
candidate family symmetries. Moreover we show that, even if we drop the constraint of an 
exact texture zero in the (1,3) position, two of the three texture zero predictions remain and 
are in good agreement with the data. Finally we prove that the identification of the phase in 
the expression for V us , eq(|2|), is true in leading order even after allowing for a matrix element in 
the (1,3) position. In |3.3| we consider the complementary possibility of an asymmetric texture. 



3 Non-zero 513 : Perturbative Analysis 

In this Section we use the notation of Hall and Rasin ||. We start with the Yukawa matrices Y 
(Y=U or D) with the assumption that the entries in the Yukawa matrices have a hierarchical 
structure, with Y33 being the largest. For the purposes of explaining the important aspects of 
the analysis it is useful to first take to be real and later consider how the analysis is modified 
by CP violating phases. The matrices Y can be diagonalized by three successive rotations in 
the (2,3), (1,3) and (1,2) sectors (denoted by S23, S13 and S12 ): 



Y n 







Y22 


(Y U 
Y 21 

\ ^31 
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^33 ) 

Y\2 1 13 

Y22 Y23 

Y32 Y33 J 



Yl 3 \ 



( 1 

s y 

V 

/ 1 





b 12 
1 




1 

6 23 







1 



b 13 



\ / 



' 23 



/ 





1 



1 





-s Y \ 


1 



V 



— s 



13 





1 





/ 1 







1 

'23 



\ 



y 

'23 



X 



' 13 



1 



6 12 



JY 
' 12 

1 







1 / 



(9) 



In terms of these angles the CKM matrix is given by 



V 



S\2 + Si 3 S 2 3 



S12 



S?3 S 23 



y _ s 13 + S 12 S 23 —s 23 — S 12 S 13 



Sl3 - sf 2 S23 
$23 + sf 2 Sl3 
1 



(10) 



where s 23 



Q D _ Q U „ 
b 23 6 23; ^13 



s 13 - s 13 and s 12 



b 12 b 12- 



It is straightforward now to see the origin of the texture zero relations eq( [I]). From eqfllDD 
we see it is sufficient to have ||: 

m \Yub\ _ | s t/ 1 an( ^ |vy. _ | s i?| ^jch j s obtained by: 



\v cb \ -1^121 m,| 



si 3 | << kf 2 S2 3 | and |si 3 | << |s? 2 s 23 |. (11) 



This condition will define how small the (1,3) element must be in the mass matrices to obtain 
the (1,3) "texture zero" prediction. Notice that since the "13" rotations are performed after 
the "23" rotations, what determines the size of the rotation si 3 is the "effective" element Y13 in 
eq. ([16|) . As eq. (|16D shows, Yi 3 depends not only on Yi 3 but also on the size of the Y 21 element, 



which is rotated in the (3,1) position by the right-handed quark rotation s'\ z . The size of the 
Y21 element can be related to the light quark masses, so that the conditions (|Tl|) also defines 
how small the right-handed rotation s'\ z must be in order to obtain the texture zero prediction. 

In addition we require 

• |s?ol = and \sP \ = which is obtained by: 

\Y n \ « |^§^| and |Y 12 | = |Y 21 |. (12) 
122 

This is the (1,1) texture zero condition together with the symmetry needed to obtain eq([TJ|2|). 
To proceed further we need to determine the mixing angles in terms of the Yukawa couplings. 
To do this we assume the off diagonal elements are small relative to the on-diagonal ones in 
each step of the diagonalisation, leading to the perturbative relation for the small mixing angles 
given by0 

Y Y23 



"2?> 



33 



Y 32 Y 22 

*33 


6 2 3 — 


Y 32 
Y 33 


Y 23 Y 22 

V 2 

r 33 


YsiYn 


b 13 — 


Y31 


Y13Y11 


V2 ' 
1 33 


Y33 


r 33 


YnYn 


JY ~ 
6 12 — 


Y21 


Y12Y11 


y2 ' 
J 22 


Y 22 


Y 2 

*22 



Y Y13 

S 13 ~ T7~ + 
1 33 

4 - ^ + ^^'12-^ + ^ (13) 

The successive rotations produce elements 

~ Y l2 Y 21 ~ Y 13 Y 31 ~ Y 23 Y 32 , 

1 11 — *11 r~ , 1\\ — J 11 77 , ^22 — 122 77 , 

122 r 3 3 Y 33 

and 

Yi 2 = Y 12 — Y^s'^, Y 21 = Y 21 — Y^s^? (15) 

Yl3 = *13 + ^31 = ^31 + ^2l4- (16) 

2 Actually this equation defines just how small the off diagonal elements need be since successive terms in 
the expansion should be well ordered. 
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From this equation one may see that the contribution of terms involving elements below the 
diagonal are suppressed by inverse powers of the heavier quark masses. For this reason they 
are only weakly constrained by the CKM mixing angles, with the only possible exception of 
the (3,2) element ||. Returning to the condition eqflTTD for the (1,3) texture zero prediction 
we see that Y13 must be small. The condition on Y31 is much weaker due to the heavy quark 
suppression. As we shall discuss the larger hierarchy in the up quark masses means that the 
down quark contribution to the mixing angles dominates. Thus, to a good approximation the 
requirement for the (1,3) texture zero prediction is 



£>13 < 



U12D23 



u- 



(17) 



22 



We are interested in what happens if this condition is not satisfied. In this case, from eqsflTOl) 



and (HID, we have 



w, 



ub 



\Vcb\ 

\Vtd\ 



\V, 



ts 




(18) 



How is this analysis affected if are complex? We discuss this in detail in the next Section 
but the implications are easy to anticipate. The sequence of rotations in eq(||) will now be 
interspersed with various diagonal rephasing matrices. This will change the above equations 
introducing phases in various terms but cannot induce any new terms in Vy. However, it is clear 
that even in this case eqs flTTD and ( [12]) are the correct conditions for yielding the predictions 
eq([l|). In turn this means that eq ([l8|) remains correct although the terms proportional to ^ 
may acquire different phases in the two equations (see below). However this does not change 
the conclusion about which texture zero prediction receives the dominant correction. 

Note that the texture zero predictions are modified in a definite way in that the prediction 



for IKJ/IK 



\Vts\ because the 



'cb\ in general has a larger percentage change than that for \Vtd\ 
mass hierarchy for the up quarks is larger than that for the down quarks while the correction 
proportional to Si 3 /s 2 3 remains the same in both cases. This is just what is needed to correct 
the disagreement we found when comparing experiment with the texture zero prediction. The 
correction term to the ratio |T4b|/|Vd,| is 

1 

— ( 

2 



1 



cu simp + -c\j 



(19) 



where c\j 

angle cb is 90° while the correction term to the ratio |Vtd|/|Vt s | is 



— l\ /— and ib is the relative phase between the two terms when the CP phase 

S23 rn " 



1 — Co cos ib H — c D 



(20) 



where c D = S -^L^. If we take $4 

S23 ' V ms l/Jool 

increase, for example, to |14& 



~ 0.04 and ib w —45° then we can easily get a 90% 
V c b\ while affecting |Vf<i|/|Vt s | by only about 10%. 
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3.1 Inclusion of phases - General parametrization 



Consider the Yukawa matrices for the case where U\\ and Du are both zero, i.e. just one 
texture zero in each. From the arguments of Kosenko and Shrock [^8| there would be then 
be a total of eight unremovable phases. We may assign the eight phases as 4>i 2 ,D , (pi 3 ,D ', 2 2 
and </> 2 3 so that all possible quantities which are invariant under re-phasing transformations, 
e.g. U12D23U22DI3 are not, in general, real - in accordance with the discussion of Kosenko 
and Shrock. Having so many phases is unnecessary as far as the physics is concerned since 



the CKM matrix elements computed from eqfllQl) in the heavy quark limit depend on only two 
independent combinations of the eight phases. Using eqs(13-16) we have 



K, 



\D 



12 



\D 



22 



\u- 



12 



\u, 



(21) 



22 



i.e. \V US \ depends on the phase 0i 



In addition we have 



V, 



ub 



til 



D 13 


. e #f 3 . 


u 12 




D 23 


D 33 




U22 




^33 



J</>23 



\D 



13 



\D 



+ 



ID 



12 



\D 



23 



^23 



33 I 



\D 



22 



\D 



(22) 



33 1 



Thus the magnitudes depend on only the combinations 2 = 

and 0i — 02 respectively. As a result so must the Wolfenstein parameters p and 77 depend only 
on the two phases 0i and 2 . Moreover, evaluating the invariant J = Im{V^&l^ s V^* V^" 6 } which 
determines the magnitude of CP violation, we find 



77 oc Im[J] 
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(23) 



For small I-D13I, the first term is the leading one and if the sub-leading corrections were neg- 
ligible, then the CP violating phase (f>cp would be just 1; the same phase which enters in 
eq(0). In this case the phase 0i is simply related to the 'standard' (i.e. PDG convention) of 
the CP-violating phase, 5by5 = 7r — 0i— f3 where f3 is the angle appearing in the unitarity 
triangle. The next leading correction is the second term, proportional to sin 02. Our proof that 
the phase 0i which drives the CP- violating phase is the same one that appears in eq(||D follows 
simply from the suppression of terms in the heavy quark limit and the equality of the (1,2) 
and (2, 1) elements. This generalizes the previous proof of this result [14], 0] which assumed an 
Hermitian form for the mass matrices with texture zeros in the (1, 1) and (1, 3), (3, 1) elements. 
Although any matrix can be made Hermitian by phase changes in general the texture zeros are 
not preserved by such transformations so this is the most general starting point. 

To summarize we have shown that it is sufficient in a parameterization of the mass matrices 
to retain two non-zero phases which we take to be 



9( 2 and 



^13 > 



i.e. 



!>f 2 and 



if 3 . Our 



analysis requires 
ratios (c.f. 0). 



90 which is the case of 'maximal' CP violation for fixed quark mass 
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3.2 Fit to the data: symmetric texture 



We now turn to a fit to the data with a non-zero entry in the (1,3) position of the down 
quark mass matrix. As we have discussed the fit to CKM matrix elements is in this case 
insensitive to the matrix elements below the diagonal. For definiteness we perform a fit making 
the assumption that the mass matrix is symmetric. This corresponds to a specific choice for 
the elements below the diagonal consistent with this "smallness" criterion. Of course, assuming 
an Hermitian rather than a symmetric form for the mass matrices does not change the quality 
of the fit. 

Leaving aside a discussion of phases for the moment, the U and D Yukawa matrices have 
the form 



(24) 



and 
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(25) 



where h t , hb are the t, b Yukawa couplings and the expansion parameters e and e will be chosen 
so that the remaining parameters a, b, c, a' , b' and d are all of 0(1). This texture is similar to 
the ansatz considered by Branco et al. . 

In fact we can extend this parameterization to include charged lepton masses by choosing 
the same matrix L (with the same parameters) for the charged leptons mass matrix as for the 
down quarks except for the usual Georgi-Jarlskog |3(J factor of —3 multiplying the (2,2) element. 
This retains mb-m T unification and also gives good predictions for the lighter generations of 
lepton. 

The small expansion parameters are determined immediately, since to leading order (all our 
discussion is to leading order) we have 



Since the above texture should apply at the unification scale, we have 



0.05 



0.15 



(26) 



(27) 



i.e. we see that e < e and suggests e = 0(e 2 ). The coefficients b and b' are also determined to 
leading order since 

(be) 2 (28) 



b'e) 2 
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giving b ~ 1.5 and V ~ 1. So far the parameters are taken real but we now introduce phases 
as discussed in Section 53 . For the case of interest with a small (1,3) element only two phases 



play a role in determining the physics, as discussed above. Here we assign a phase to each of 
U and D by taking 

b' -> 6' e** c -> c (29) 

Having chosen to attach a phase to the (1,2) element of U, as discussed in section [O] the 
one phase chosen in D should be attached to a different element of D. We then expand all 
the elements of the rotation matrices which diagonalize U and D in terms of e and e, retaining 
the leading terms only. The CKM elements are expressed in terms of these rotation matrix 
elements (see |§ for example) which then allow leading order expressions for the V^j KM in 
terms of our expansion parameters. 

At this leading order, we have 



V us = be-b'ee i4 > \V US \ = A = J— (1 + 0{e/e)) (30) 

V m s 

where A is the Wolfenstein parameter. The Wolfenstein parameter A fixes the value of \V cb \ 
and in our expansion 

\V cb \ = AX 2 = at + 0(et) =^ A = £ . (31) 

We are neglecting here the up quark contribution to \V cb \, which is justified if the (2,3) elements 
in the up quark matrix are indeed of order e 2 as suggested by eq. (f23j) . From a phenomeno logical 
point of view, however, a larger size for those elements (e.g. of order e) is also allowed and 
could lead to non negligible up quark contributions to V cb . 

Since the textures that we are discussing apply at the unification scale rather than at low 
energies, we must use values for mass ratios and CKM parameters appropriate to that scale. 
For this it is convenient to introduce the parameter x — (Mx / Mz)~ h * ^ 16tt2 ^ ~ 0.7 and then 



x, 



A(M Z ) 
{m s /m b ){M x ) 
(m s /m b )(M z ) 

(m c /mt)(Mx) 
(m c /m t ){M z ) 



X 
X 

3 



r (32) 



Thus at the unification scale we have A ~ 0.58 or a ~ 1.3 As discussed above, up to corrections 
suppressed by inverse powers of the third generation masses, the phase determines the sign 
and magnitude of the CP-violating CKM phase. In the SM context, the observed CP violation 
requires a near maximal phase, <fi ~ 90°, so 

V ub = cte^-iab'et 

V td = -ceV^ + abe 3 (33) 
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which imply that, to leading order the Wolfenstein parameters which govern the size of V ub and 
V td are given by 




e , b'e . , 

H — - c cos ip — c sm ip 

ab ab 2 



ab 



csini/j 



ab 2 



c cos ib 



(34) 



For c small, the phase <fi ~ +90° fixes the correct sign of the first (dominant) term in the 
expression for 77 and maximizes CP violation for fixed quark mass ratios. In passing we note 
that, to this order, the entire list of quark masses and CKM matrix elements do not involve a' 
or d and so could take on any value of 0(1) without affecting the physics. We can express the 
perturbation to the canonical values for |\4b/K;fe| and |Vt<f/Vts| given by eq([|) 
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(35) 



The effect of filling in the (1,3) texture zero of D can be more dramatic for |V^b/V^,| than for 
|Vfd/Vta| since the correction to the latter is suppressed by e/e ~ \Jm u /m c / ' \Jmd/m s relative 
to the correction to the former. We can get the desired phenomenological result of moving 
iKih/K&l up towards the measured value around 0.09 while not unduly perturbing the value of 
\V t d/V ts \ given by eq(0). 



Parameter 

Q 

m u /m d 
m c /m s 



Value 
22.7 ±0.8 
0.533 ±0.043 
9.5 ± 1.7 




Table 3: Values for the quark masses ratios and the parameter Q (defined in eq([36|)) used for 
the texture zero fits. 



Using the expansions of eqs (|24]j25|) for the U and D matrices in terms of e and e we carry 
out a fit using information on the measured estimates for ratios of quark masses and CKM 
matrix elements. The expansion parameters e, e and the 0(1) coefficients a, b, b' and c are 
determined via the expressions ( |2^J2^J30|J3T|J34]J35| ) . We follow BHR in using the combination 
Q given by 

m s /m d 



Q 



(36) 



y/l - {m u /m d ) 2 

which is determined accurately from chiral perturbation theory. Additionally we can use the 
ratios of the masses m u /m d , m c /m s . 
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Figure 3: Fit A of Section 3.2 (symmetric texture) to the measurements of |Kt&/Vd>|, Am# s , Am^, 
| | and sin 2(3. 

The resulting fit (Fit A) yields the values 

e = 0.15 ±0.01 |6| = 1.5±0.1 a = 1.31 ±0.14 
\c\ = 2.7 ±0.10 ^ = -24° ±3° (37) 

where the values are given at the unification scale. 

Demanding a texture zero in the (1,3), (3,1) elements for a symmetric D matrix leads, in 
particular, to too small a value for p and there is a marked improvement in the overall description 
of the data when this zero is filled in, as can be seen by comparing Fig. ^| with Fig. |2|. 

That c ~ 3 means that the order of the (1,3) term is ambiguous and could be either 0(e 4 ) 
or 0(e 3 ). Q While the texture zero in D 13 , D 3 i does lead to the rather attractive result of eq(|l]) 
the current experimental data, within the context of the Standard Model, favour perturbative 
corrections as suggested by eqs(^). 

We noted earlier (Section |2.1|) that there is a potential disagreement in the SM fit to |e^|, 
Amg s , Am^ and the recent measurement of sin 2/3 = 0.41 ± 0.170. To quantify the impli- 
cations of this, we perform a second version of Fit A (still using the texture of eqs (p4]j25|) ), 
dropping the constraints of |e^|, Amg s and Am^. This is shown in Fig. f|. The effect is to 
increase slightly the value of the parmeter c to 3.31 ± 0.10 and the value of ip to 6° ± 3°. 

In Fig. ^| we show the resulting probability distributions for the quantities y^, \ck\i Am Bs 
and sin 2/3 compared to the corresponding experimental distributions. 

3 That d is quite undetermined means that a texture zero in U13, U31 is not ruled out. 
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Figure 4: Fit A of section 3.2 (symmetric texture) to the measurements of IV^/V^I and sin 2/3. 

There is actually a further solution (Fit B) to the above equations where the phase ~ —90° 
rather than +90° which corresponds to to V — > —f in eqs( |3"4] , |3~5|) . In this case the first term in 
the expression for rj is no longer the dominant one and a larger value of c is needed to obtain 
a positive value for rj. Only the parameters c and ip change from Fit A, and for Fit B we find 
c = 8.45 ±0.33, ip = —58° ±5°. The result of Fit B is shown in Fig. |B| and comparing with Fig. [5] 
we see no difference in the quality of fits A and B. This solution has the (1,3) matrix element 
of the same order (0(e 3 )) as the (1,2) matrix element. In this case Di 3 /Di 2 ~ D 23 / D 22 ~ 1 ; 
suggestive of a non-Abelian family structure. 

3.3 Fit to the data: asymmetric texture 

We now consider a fit to the data in which the texture zero relation = ^/m u /m c 

is modified by higher order corrections in the perturbative diagonalization of the down mass 
matrix. As in the previous section, the correction to that relation comes from a small but 
non negligible rotation sf 3 induced by a non zero element -D13. Here, however, we assume that 
-D13 is mainly induced by the rotation s' 23 used to diagonalize the 23 sector of the down quark 
mass matrix, the initial value _D 13 being negligible (see eq. flT6|)). This situation is therefore 
complementary to the one considered in the previous subsection, where Di 3 was mainly given 
by the original entry Di 3 and the contribution to Di 3 proportional to s' 23 was assumed to be 
negligible. 

Let us first of all estimate how large s' 23 should be in order to give a significant contribution 
t° \V u b/Vcb\- We start from eqs. (|T8|), that assume only |F 12 | = \Y 2 i\, Y n = 0. The size of the 
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Figure 5: One dimensional probabilities for the physical observables IK^I/IV^I, Arri£j d and sin 2/3. 
The probabilities in red (dark) correspond to the experimental constraints and the probabilities in 
blue (light) correspond to the predictions of the texture with 13 and 31 entries different from zero 
with c = (2.7 ± 0.10) and ip = (-24 ± 3)°. 



correction is determined by s^, which can be written as 

Y Yi 3 \m{m 2 lY fQQ , 
*13= xT~\ — — s 23 > (38) 

-i 33 V m2 777,3 

where m» is the mass of the quark of the i-th family in the sector Y = U,D and we have used 
Y13 = F^s^j Y12/Y22 ~ <Jm\jm 2 and ^22/^33 — rn 2 /m 3 . From eq. one can see that the 



contribution of the up quark rotation to s 13 = sf^ — sf 3 is negligible. The mass ratios in 
eq. ( j38|) are in fact much smaller in the up sector than in the down sector. Moreover, the 
factor s' 2 \ in (j3"8"|) appears in the product s^s^, which contributes to the ratio m 2 /m 3 . Barring 
cancellations, we therefore have s' 2 \s 23 <m 2 /m 3 , a constraint stronger in the up sector due 
again to m c /m t <C m s /mb. Hence we can safely neglect sf 3 and write 

S13 S13 [m^m s /m b 



S23 S23 V m s \V C 



cb\ 



=)' D 
'23 



(39) 



Since at the electroweak scale we have ^Jmd/m s ~ 0.22 and (m s /mi)/\ V c b\ ~ 0.6, in order to 
get a correction s i3 /s 23 ~ 0.03 (which gives a good agreement with data in absence of phases) 
it is sufficient to have s'23 ~ 0.2 — 0.3, or \D 32 \ ~ (0.2 — 0.3)|D 33 |. Larger values are also in 
principle allowed depending on the phases in eq. (|T9"|). Notice that in first approximation the 
ratio Z} 32 /.D 33 at the unification scale is the same as at the electroweak scale. 

Quark mass textures with \D 32 /D 3S \ = 0(1) have been considered in the literature Q in 
connection with a large leptonic mixing angle originating from the charged lepton mass matrix. 
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Figure 6: Fit B of Section 3.2 (symmetric texture) to the measurements of |Fu&/V^,|, Atob s , ArriB d , 
| ex | and sin 2(3. 

In SU(5) unified models, the ratio I-D32/ -D33I associated with the right-handed quark rotation 
s'i 3 corresponds in fact to the ratio of charged lepton mass matrix elements \E2s/E^\ associated 
with a left-handed charged lepton rotation that mixes /1 and r neutrinos. On one hand, this link 
can be considered as a motivation for studying asymmetric textures as a solution of the V u b/V c b 
problem. From the opposite point of view, we can say that in the context of unified models, the 
asymmetric solution of the V u b/V c b problem, pointing at a largish value of I.D32/-D33I, indicates 
that the large mixing angle responsible for the atmospheric neutrino anomaly comes, or at least 
receives a significant contribution, from the diagonalization of the charged lepton mass matrix. 

The precise relation between the right-handed quark and left-handed lepton rotations de- 
pends on possible Clebsh coefficients relating transposed down quark and charged lepton matrix 
elements. The presence of non-trivial coefficients enhancing some charged lepton matrix element 
is indeed suggested by the empirical relation m^/m T ~ 3m s /m b . Since a value of I.D32/-D33I 
around 1/3 is preferred by a fit of data, a near to maximal lepton mixing can be obtained if 
the Georgi-Jarskog factor 3 sits in the E23 entry. 

Notice that a sizeable |-D 32 | ~ | -D33 1 / 3 indicates an asymmetry I-D32I | -D23 1 • The ratio 
I-D23/-D33I is in fact expected to be of order \V c b\ <C 1/3, barring cancellations between up and 
down quark contributions to V c b- Such an asymmetry can be easily obtained in the context 



of Abelian and non-Abelian models. In Section [4.2| we will describe an explicit example of 
non-Abelian family symmetry leading to the asymmetry in the 23 sector while preserving the 
relation |y 12 | = |Y" 2 i| m the 12 sector and the texture zeros Yi 3 ~ y 31 ~ Y\\ ~ 0. This texture 
allows us to isolate and study the corrections to the texture zero relations we are considering 
in this section. 
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To examine the implications of this scenario we modify the parameterisation of the Yukawa 
matrices given in eqs(|24|, |25| ) by restoring the exact texture zeros in the (1,3) and (3,1) elements 
and allow \D 32 \ 3> | 1 • 

We consider the following parameterisations for the mass matrices for the up and down 
quarks whose absolute values can be parametrized as 
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(40) 



(41) 



with < e < 1. In this parameterisation of D we see, comparing with eq ([25D , we have 
changed the (3,2) element to be 0(1). The remaining elements are of the same order, e in 
eq(|41|) being 0(e 2 ) and e' is 0(e 3 ). The parameterisation of U has the same form as eq(f2~4|) 
with the exception of the (2,3) and (3,2) elements which are now of 0(e) and not 0(e 2 ). We 
have chosen this form to relate to a promising texture model discussed in Section [4~2"| . However 
these elements are poorly determined by the data and it is possible to obtain solutions where 
a, b are O (e), corresponding to the original symmetric parametrisation of eq(|24]). 

In order to obtain the precise form of the corrected texture zero relations, we first write the 
general expression for \V US \, \V u b/V cb \, \V t d/Vt s \ in terms of the rotations defined by eq. (01): 
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(42) 



where t® 2 and tf 2 are the tangent of the 12 rotation of left-handed down and up quark re- 
spectively. Written as above in terms of the tangent of the angles (and cosines c 



D 
12) 



-127) 



the 

expressions are exact up to O (A 4 ) corrections, A being one of the Wolfenstein parameters. The 
phases 0i, 



i>2 were discussed in Section 3.1 



In order to obtain the relation between mixing 
angles and quark mass ratios generalizing the texture zero relations, one has to express the 
angles in eqs. 



2|) in terms of quark masses. In the case of textures ( |4l] , f40"D , we obtain 
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(43) 



where t = ID32/-D33I as in eq. (f£T|). The parameter t represents the tangent ti£ of the 23 
rotation on the right-handed down quarks in eq. <H). Since we are considering the possibility 
of a sizeable t, we do not approximate the cosine of that angle, c = l/y/l + 1 2 , with 1. As a 
consequence we have, at leading order in A 2 = O (md/m s ), 



m q V c m q 



This is because the diagonalization of the 23 sector in the down sector not only induces an 
effective (1,3) element, but also generates a slight asymmetry in the 12 sector: \Di 2 \ 7^ I-D21I 
despite \Di 2 \ 7^ | Z? 2 i | - The expression for tf 2 above also includes a next to leading correction 
in A 2 = O (md/rris). 

Let us now turn to the numerical determination of the parameters entering the expressions 
for U and D. Three of those parameters, t and the two phases <f>\, <p 2 , can be determined 
independently of the values of the others in terms of \V US \ and p, fj (as obtained from the SM 
fit). This is possible since, for given values of the quark masses, eqs. ( f42| , [43"D relate t, 0i, <p 2 
to |V^ S | and IV^/V^I, | V^^/ V* s | , and therefore to \V US \ and p, fj. More precisely, to get the best 
values of t, <pi, 2 , we use the following procedure. First we calculate (K^/V^l, |Vtd/Vt s | in terms 
of p, fj. Then, for any given value of the ratio t = I.D32/-D33I and of the phasef] <f> 2 we recover 
md/m s and m u /m c from eqs. (^,^) (the phase 0i is obtained, up to discrete ambiguities, from 
the relation involving |V^ S |). Finally, we calculate Q, m u /md in terms of md/m s , m u /m c . We 
can at this point perform a fit of Q, m u /md in terms of p, fj for any given value oft and <p 2 . The 
quantities m c /m s , m s /mb, \V c b\, A involved in the relation between Q, m u /md and p, 77 are also 
included in the fit. As mentioned before, using Q, m u /md instead of md/m s , m u /m c improves 
considerably the quality of the fit, especially for small t (despite it requiring the inclusion of 
the ratio m c /m s ) [n]. We then obtain t = 0.3, <p 2 ~ — 0.2tc, cos0i ~ 0.1. Notice that 0i again 



turns out to be almost maximal as was the case with symmetric textures. 

The determination of t, 1; 2 described above depends on the additional parameters 
in ([K],[n]) only through the quark masses and \V c b\. As a consequence, the fit of all data 
in terms of all parameters decouples into a fit of \V US \, p, f] in terms of t, <fii, <p 2 (which makes 
use of the experimental values of the quark masses and |K&|) an d a fit of quark masses and \ V c b\ 
in terms of the additional parameters. The first fit has been described above and is independent 



of the specific form of the textures (||,40) and of the values of the additional parameters one 
uses to account for the quark masses and \V c b\. For example, independently of whether the up 
quark matrix is fully symmetric or not, the parameters a and b are of 0(1) or smaller, a, (3 are 
of 0(1) or vanishing, provided that the values of the quark masses and of \V c b\ can be accounted 
for. All the dependence on the specific form of ( fTO||IT| ) is therefore confined to the fit of quark 
masses in terms of the parameters e, e', a, b, c, a, (3. Here we consider in more detail such a fit 



in the case motivated by the flavour model described in Section |4.2| in which a — (3 — 0. 



To be precise, in order to include all sign ambiguities in a single phase we actually fix the value of the phase 
b' 2 defined by cos0' 2 = cos 02, sin^J, = sin 02sign(sin In the following, 4>2 should be read <p' 2 . 
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Let us start from the quark masses. Besides m t , m&, trivially accounted for by h t} h^, we 
have to account for m s /m{,, m c /m t , m^ms/m 2 ,, m u m c /m 2 . We work at the m t scale. In order 
to account for m^/m^ ~ t e we need e ~ 0.08. From m c /m t ~ abe 2 we then have ab ~ 0.6. Since 
m,dm s /ml = e /2 /(l + t 2 ) 3 / 2 , we can obtain e' = 0.006. The ratio m u m c /m 2 = (cee') 2 then gives 
c ~ 0.5. Notice that all parameters that are supposed to be of 0(1) indeed are. In particular, 
we have \U 2 3/h t \ ~ \U^/h t | ~ \D 2 2,/h h \, which justifies the use of the same parameter (e) in 
both the U and D matrices. In the model of Section |4.2| the same parameter indeed appears in 
those entries because the same vev generates them. Finally, we have to account for the value of 
\V c b\, which has a contribution from the down sector, e, and one from the up sector, be. Both 
contributions are of the right order of magnitude, so it is clear that \V c b\ can be obtained for 
an appropriate choice of the 0(1) coefficient b. The precise value of b depends on the relative 
phase between the two contributions. For a = (3 = 0, the phases of the SM quark multiplets 
can be redefined in such a way that U/ h t and D/hf, differ from their absolute values in ( pEl|J£0| ) 
only by a phase e*^ multiplying t in D32 and a phase e 1 ^ multiplying b in f/23. The phases (f)\ 
and 02 are then given by 



e i<t>i = e W-4>) (44) 
e (be^ - l) 

„*h - >. . (45) 



The relative phase between the two contributions to turns out to be e 1 ^ and can be de- 
termined from the equations above. We are actually interested to the value of b which simply 
follows from eq. (flop: b = \e l ^ 2 \V c b\/ e + 1| ~ 1.4. We then also have a ~ 0.4. Notice that a 
and b are both of 0(1) but b > a. This slight asymmetry is reduced by the running up to the 
unification scale. 

The results for this fit with the asymmetric texture, D32 3> -D23 are very close to those for 
the symmetric texture of 3.2 and the resulting contour plot in the p - fj plane is essentially 
identical to that in fig. [| 



4 Implications for a family symmetry 

Of course the underlying motivation for studying the detailed structure of the quark and lepton 
mass matrices is that they may lead to an insight about the structure beyond the Standard 
Model. Here we briefly comment on the implications of our analysis for such structure con- 
centrating on the possibility there is an extension of the symmetries of the Standard Model to 
include a family symmetry. 



4.1 Symmetric case 

It turns out to be remarkably easy to construct a model leading to the mass matrices in 
eqs. ([24],[25|) through the introduction of an Abelian gauge symmetry, U(l) (such additional 



22 





Q t u\ d\ Li e£ if H 2 H x 


U{l) FD 


OLi OL{ OLi di CLi di — 2«x — 2«i 



Table 4: U(1)fd symmetries. 

symmetries abound in string theories). The most general charge assignment of the Standard 
Model states is given in Table This follows since the need to preserve SU(2)l invariance re- 
quires (left-handed) up and down quarks (leptons) to have the same charge. This together with 
the requirement of symmetric matrices then requires that all quarks (leptons) of the same i-th 
generation transform with the same charge ajj(aj). If the light Higgs, H 2 , Hi, responsible for the 
up and down quark masses respectively have U(l) charge so that only the (3,3) renormalisable 
Yukawa coupling to H2, Hi is allowed, only the (3,3) element of the associated mass matrix 
will be non-zero as desired. The remaining entries are generated when the U(l) symmetry is 
broken. A particularly interesting example may be constructed in a supersymmetric extension 
of the Standard Model []3l]. We assume this breaking is spontaneous via Standard Model singlet 



fields, 9, 9, with U(1)fd charge -1, +1 respectively, which acquire vacuum expectation values 
(vevs), < 9 >, < 9 >, along a "D-flat" direction. After this breaking all entries in the mass 
matrix become non-zero. For example, the (3,2) entry in the up quark mass matrix appears at 
0(e' Q2 ~ ai ') because U(l) charge conservation allows only a coupling c c tH 2 (9 / M 2 ) a2 ~ ai , a 2 > cti 
or c c tH 2 (9 1 M 2 ) ai ~ a2 , oil > ci 2 and we have defined e = (< 9 > /M 2 ) where M 2 is the unification 
mass scale which governs the higher dimension operators. As discussed in reference |3l|] one 
may expect a different scale, Mi for the down quark mass matrices (it corresponds to mixing 
in the H 2 , Hi sector with M 2 , Mi the masses of heavy H 2 , Hi fields). Thus we arrive at mass 
matrices of the form 



Mu 
m t 



( /HiPn4 2+6a| h 2 pi 2 e[ 3al hi 3 p n el^ \ 

h 2 ip 2 ie[ 3 ' h 22 p 22 e 2 h^p^e 1 
V />3iP3iel 1+3a| h 32 p 32 e l h 33 J 



(46) 



I fc n <7nel 2+6a| ki 2 ai 2 e b W ki 3 cn 3 e a ^ \ 

^2iO"2iefe |3a| k 22 a 22 e 2 k 23 a 23 ^ (47) 
y /c 31 cr 31 e a l 1+3a l k 32 a 32 e l k 33 j 

where e = (^)' Q2 " ai1 , e = (^) |Q2_ai1 , and a = (2a x - a 2 - a 3 )/3{a 2 - at x ). For -3a > 1 
e a = e& = e and e a = e b = e. In this case it is easy to check that there are no texture zeros 
because all matrix elements contribute at leading order to the masses and mixing angles. For 
1 > —3a > 0, e a , e a change and are given by e a = (^f^)'" 2- " 1 ', e a = (ivf^)'" 2- " 1 '- m this case 
texture zeros in the (1,1) and (1,3) positions automatically appear for small < 9 >. However 
the (1,2) matrix element is too large (cf Table 0). For a > however e a ,b = (^f^)' 02- " 1 '; 
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e a,b = (ivlf)' 02 ai 'i texture zeros in the (1,1) and (1,3) positions persist, and the (1,2) 
matrix element can be of the correct magnitude. 

Note that the family symmetry does not make the small elements exactly zero so it predicts 
only approximate texture zeros. Indeed, fixing the parameter a = 1 to obtain the measured 
magnitude of the (1,2) matrix element one finds that the (1,1) elements occurs at 0(e s ) and 
0(e 8 ) for the up and down mass matrices respectively. This is so small that eq(fj|) is valid to a 
high degree of accuracy if pu = pi\. On the other hand the (1.3) matrix element is predicted 
to occur at 0(e 4 ), 0(e 4 ) for the up and down matrices respectively and, as discussed above, a 
term of this order (with a coefficient 2) is sufficient to correct the prediction for ^ following 
from the assumption of an exact texture zero. 

However the best fit prefers the (1,3) element to occur at 0(e 3 ), 0(e 3 ), i.e. close to the 
(1,2) matrix elements. Moreover the measured value of requires the (2,2) and (2,3) matrix 
elements of the down quark mass matrices should be of the same magnitude, of 0(e 2 ). This is 
in contradiction to the predictions of the Abelian family symmetry, unless one appeals to the 
unknown coefficients of 0(1). The most plausible way to get such relations for matrix elements 
involving different family members is to invoke a non- Abelian family symmetry H. 



4.2 Asymmetric case 

We now describe a supersymmetric non-Abelian model based on a U(2) family symmetry acting 
on the two lighter families [f27| and on the unified gauge group SU(5). This model is a varia- 
tion which leads to asymmetric textures discussed in Section |373] with Y" 2 2 ~ 0. The lighter 
families ip a , a = 1, 2 (ip = T, F, where T and F are respectively the 10 and 5 representations of 
SU(5)) transform as ip a — > U a b^b under U GU(2), whereas the third family and the Higgs fields 
Hi, H 2 are invariant. Such a symmetry is approximately realized in nature. In fact, in the U(2) 
symmetric limit the lighter fermion families are forced to be massless and in supersymmetric 
models their scalar partners are forced to be degenerate. The symmetry is broken by two SM 
singlet scalars, an antidoublet <p a transforming with U T ~ l and an antisymmetric tensor A ab 
transforming with U T ~ l <8> U T ~ l under U(2). In an appropriate basis in the flavour space, the 
corresponding vevs can be written in the form 



W = : W = .. . (48) 





where V, v > 0. The correct hierarchy and mixing between the two lighter families is obtained 
if v/V = 0(\V US \). The U(2) breaking is communicated to the light fermions by an heavy 
U(2) anti-doublet x a through a Froggatt-Nielsen mechanism. Under the gauge group, each x a , 
a = 1, 2, transforms as a full fermion family, which allows a mixing with the light fermions. The 
heavy mass term Mx a Xa for the fields x a a l so involves of course a doublet Xa with conjugated 
tranformations under the SM and U(2) group. As for the size of the mass term, one simple 
possibility is that the scale M is above the SU(5) breaking scale, M > Mgut- A small ratio 
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V/M is then generated if the U(2) breaking takes place at the SU(5) breaking scale, V ~ 
Mqut- Possible SU(5) breaking corrections to the heavy mass M will also be correspondingly 
smaller. The small ratio V/M determines the small Yukawa couplings accounting for the second 
generation masses and mixings, forbidden in the U(2) symmetric limit. In particular, the correct 
order of magnitude for the mixing of the two heavier families is obtained if V/M ~ O (\V c b\). 

Besides x a ■, Xa, representing the minimal choice for the messenger sector, the physics at the 
GUT scale can involve additional heavy fields. For example, we have mentioned in the previous 
Section the possibility of a mixing in Higgs sector involving two heavy Higgs fields H[, H 2 , 
singlets in this case under the U(2) symmetry as Hi, H 2 . Such a mixing can be used to account 
for the hierarchy mj <Cm t . We therefore include H[, H' 2 in the model. Since they are allowed 
to interact both with the light families and the U(2) breaking sector, the U(2) singlets H[, H' 2 
can also mediate U(2) breaking. Notice that the scale M' at which this singlet-mediation takes 
place is a priori independent of the scale M associated to the doublet-mediation. For example, 
if the mass of the heavy Higgses is set by SU(5) breaking we will have M' <C M. 

At this point one can write the most general renormalizable superpotential involving the 
light fermions (ip a , ips), the Higgs fields (Hi, H 2 ), the U(2) breaking fields (<p a , A ab ) and the 
doublet (x a , Xa) arid singlet (H[, H 2 ) messengers. Once U(2) is broken, a mixing between the 
previously massless fermions and the heavy messengers is generated. The new light fermions 
can then be easily identified by diagonalizing the heavy mass matrix. This leads to the following 
textures for the up and down quark mass matrices: 



D/h b 



U/ht 



J 



( e' \ 
-e' e 
t 1 

cee' \ 
-cee' fee 
y ae 1 J 



(49) 



(50) 



where e = O (V/M), e' — O (v/M), t = O (V/M') and all other coefficients arise from couplings 
of order one. One then obtains the textures ( |4"T| , |4"0"| ) for the absolute values of the mass matrices 
with a ~ (3 ~ 00. In particular, the relation \Di 2 \ = \D 21 \ follows from the symmetry properties 
of the U(2) representations^]. 

A few comments are in order. Since we expect \D 23 /D 33 \ = 0(\V c b\), the texture zero 
in the (2,2) position requires t = \D 32 /D 33 \ ^> \V C0 \ in order to account for the value of 
\D 32 /D 33 ■ D 23 /D 33 \ ~ m s /mb = O (\V c b\). We therefore expect a sizeable t>6, which leads to 
a non negligible correction to From the model building point of view, the asymmetry 

I-D32I 3> I-D23I corresponding to t ^> e can be simply accounted for by a relatively light singlet 
messenger scale M' <C M (which is analogous to the M\ -C M 2 assumption of the Abelian 

5 Unlike in eqs. 

here the O (1) coefficients can be complex. 
6 A non negligible correction to \D\%\ = \D%i\ can arise if t = 0(1) from the diagonalization of the kinetic 
term. 
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case). In fact, if this is the case, the leading contribution to D 32 comes from the exchange of 
the U(2) singlets H[, H' 2 at the scale M' . As a consequence, | -0 32 j turns out to be larger than 
| Z?23 1 , which is generated by the exchange of the U(2) doublets x a ■, Xa at the higher scale M. 
Moreover, H[, H' 2 transform as 5 and 5 of SU(5). Therefore, the singlet exchange at the lower 
scale M' does not contribute neither to the (2,3) nor to the (3,2) element in the up quark mass 
matrix, so that both {7 2 3 and U32 are of order e. The larger hierarchy m c /m t <C m s /mb follows. 
As for the further suppression of m u m c /mf with respect to m d m s /ml, that is automatically 
acheived if A ab is a SU(5) singlet. The operator A ab T a T b H vanishes in fact in this case due 
to the antisymmetry of A ab . This is a generic appealing feature of U(2) models. In order 
to generate a non- vanishing Uu entry, SU(5) breaking effects must be included either in the 
messenger masses or through higher dimension operators, thus giving the extra e in Uu- 



5 Conclusions 

The presence of texture zeros in the quark Yukawa matrices can constrain quite tightly the 
detailed features of the CKM unitarity triangle. Recent data has shown that it is no longer 
viable for s 13 to be zero and are, as a result, inconsistent with the most promising texture zero 
structure. This result seems quite reliable, following both from the improved bound on Atubs 
as well as the improved value of |14&/Vd>|- 

Theories invoking family symmetries beyond those of the standard model can lead to a 
hierarchal structure for U and D in which the elements appear in the form e k , where e is a 
small parameter. Motivated by such an expansion, we have explored a perturbative approach 
in which the rotation S13 is small but non-zero. One way to do this is by allowing small entries 
to replace some of the texture zeros. We have investigated a common symmetric form for U 
and D where, in particular, the (1,3), (3,1) element is no longer zero. We have also investigated 
an alternative possibility for generating S13 by allowing an asymmetric form for the (2,3), (3,2) 
mass matrix elements. Both cases lead to a desir able ph enomenological result whereby the 
perturbation of the ratio |V^&/Vd,| from the value ^m u /m c is larger than that of |V^/V ts | from 

^Jm d /m s . 

On the theoretical side, we would hope that pinning down the allowed structures for the 
Yukawa matrices will provide clues to the nature of an underlying family symmetry. Our 
analysis shows that a perturbative expansion in terms a small parameter is quite successful 
and therefore supports the idea that the family symmetry is spontaneously broken at the high 
energy scale. This has a concrete realization in the Froggatt Nielsen mechanism where light and 
heavy states are mixed via an extension of the 'see-saw' mechanism. If we require the symmetric 
form of the mass matrices it is necessary to have non-vanishing (1,3), (3,1) matrix elements. 
This is an interesting result because it follows from specific Abelian (and non-Abelian) family 
symmetries. Similarly the asymmetric solution can also be obtained from non-Abelian family 
symmetries. Improvements in the measurements of the quark masses and CKM mixing angles 
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and in particular on sin 2(3 ( as well as on sin 2a) should help in distinguishing between these 
candidate symmetries and possibly lead to a viable theory of fermion mass generation. 
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